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Resumen 


La operación y seguridad de presas han sido de gran importancia en 
México durante muchos años. Las políticas de operación en presas 
son ampliamente estudiadas en la hidráulica moderna. Los eventos 
extremos que afectan al país provocan que dichas políticas se revisen 
año con año. Este trabajo presenta una aplicación de la técnica 
multivariada de análisis canónico para correlacionar los valores de 
intensidad de precipitación durante un evento huracanado, calculados 
con la técnica del hidroestimador, con la variación de los volúmenes 
dentro de un embalse, sin utilizar características paramétricas de las 
cuencas de aportación. Se presenta el caso del huracán Carlotta 
(junio 2012), obteniendo la estimación del volumen de los 
almacenamientos en las presas Ing. Carlos Ramírez Ulloa, Valerio 
Trujano, Tepuxtepec e Infiernillo. Los resultados muestran una 
correlación superior a 0.90 en la relación existente de las 
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intensidades de precipitación huracanada y la evolución de los niveles 
en presas. Estos resultados permiten considerar al análisis canónico 
como una alternativa para detallar las políticas de operación de 
presas en México. 


Palabras clave: correlación canónica, variables  canónicas, 
eigenvalor, huracanes, presa, volumen de almacenamiento 


Abstract 


The operation and the security of dams have been of great 
importance in Mexico for many years. The policies of operation in 
dams are widely studied in modern hydraulics. The extreme events 
that affect our country cause these policies to be reviewed year after 
year. The present work presents an application of the multivariate 
technique of Canonical Analysis to correlate the values of intensity of 
precipitation during a hurricane event, calculated with the technique 
of the hydrostatimator, with the variation of the volumes inside a 
reservoir, without using parametric characteristics of the contribution 
basins. The case of Hurricane Carlotta (June 2012) is presented, 
obtaining the estimation of the volume of storage in the dams: Ing, 
Carlos Ramírez Ulloa, Valerio Trujano, Tepuxtepec, and Infiernillo. 
The results show a correlation higher than 0.90 in the existing 
relationship of the hurricane rainfall intensities and the evolution of 
levels in dams. These results allow us to consider the Canonical 
Analysis, as an alternative to detail the policies of operation of dams 
in Mexico. 


Keywords: Canonical correlation analysis (CCA), canonical variables, 
eigenvalue, hurricane season, dam, storage volume 
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Introducción 


El análisis de correlación canónico (canonical correlation analysis, 
CCA), es una técnica utilizada en casi todas las disciplinas 
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relacionadas con los recursos hídricos. En la actualidad, el CCA se 
utiliza para describir la relación que existe entre variables 
climatológicas, fisiográficas y ambientales que intervienen en el 
proceso lluvia-escurrimiento en el entorno regional de una cuenca 
(Barnston € Ropelewski, 1992; Roundy, Yuan, Schaake, € Wood, 
2015). Este complejo mecanismo de respuesta de una cuenca no es 
fácil de estudiar. Para conseguir un íntegro conocimiento de los 
procesos del ciclo hidrológico es frecuente llevar a cabo estudios de 
regionalización hidrológica, lo que permite cuantificar la correlación y 
homogeneidad entre las variables mencionadas (Kamali-Nezhad et 
al., 2010; Archfield, Pugliese, Castellarin, Skoslashien, € Kiang, 
2013). El CCA, al igual que otras técnicas multivariadas empleadas en 
la hidrología, como el análisis de componentes principales (empirical 
orthogonal — function, EOP); análisis de correspondencia 
(correspondence analysis, CA), y el análisis factorial de 
correspondencia (correspondence factorial analysis, AFC), permiten 
conocer no sólo la correlación y homogeneidad sino también la 
estructura de las variables contenidas en un modelo hidrológico 
dentro de un espacio invariante (Yu, Chu, € Schroeder, 1997; Shu 8 
Ouarda, 2007). Estas técnicas estadísticas son de gran utilidad en la 
hidrología moderna y se basan en el cálculo numérico de un conjunto 
de variables, por ejemplo, caudales y niveles, representadas en un 
modelo matemático a partir de una combinación lineal de variables 
independientes, como lámina de lluvia, intensidad de precipitación, 
temperaturas, etcétera (Steinhorst 8 Williams, 1985). La ventaja 
principal de esta formulación es que no hay necesidad de establecer a 
priori qué combinación de variables independientes se relacionan con 
los valores de entrada (Sharif, Burn, 8 Hofbauer, 2013). La mejor 
combinación lineal de la mejor correlación se encuentra durante el 
desarrollo de la metodología del CCA. Asimismo, el análisis define, de 
entre todas las variables que miden el rendimiento del sistema, cuál 
es la mejor forma de combinarlas linealmente (Fiering, 1982). 


Los modelos de lluvia-escurrimiento han sido ampliamente utilizados 
como herramientas para la predicción de caudales torrenciales en 
cuencas urbanas y cuencas de gran capacidad de drenaje (Licznar, 
De Michele, 8 Adamowski, 2015). Recientemente, la evaluación del 
impacto ambiental que sufre una cuenca hidrológica, producto de un 
posible cambio climático, es también estudiado, evaluando las 
respuestas del proceso lluvia-escurrimiento (Marzban, Sandgathe, 4 
Doyle, 2014). Hoy día, la respuesta de una cuenca a la entrada de 
precipitación es de gran importancia y se utiliza en la ingeniería para 
el diseño de obras hidráulicas, como presas, por ejemplo, para el 
abastecimiento o control de avenidas. Por lo general se sigue un 
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procedimiento de dos etapas: parametrización del modelo utilizando 
un conjunto de observaciones meteorológicas históricas y datos de 
escurrimiento a través de la calibración del modelo; así como la 
implementación de un modelo calibrado para predicción de la entrada 
de avenidas en el embalse (Hundecha, Ouarda, 8. Bárdossy, 2008). 
En años recientes, en México, este procedimiento de estimar 
avenidas de ingreso a un almacenamiento ha sido estudiado 
ampliamente; incluso cálculos con un tránsito inverso en el embalse 
han servido para verificar eventos de diseño. Derivado de estos 
estudios se han puesto en marcha algoritmos para determinar la 
función de respuesta en una cuenca, y con ayuda de la 
parametrización de modelos lluvia-escurrimientos proponer políticas 
de operación en grandes presas en el país (Morales-Velázquez, 
Aparicio, € Valdés, 2014). 


En México existen más de 5 000 presas y bordos. Si se concentra el 
volumen de agua en las 180 presas más grandes del país, se tendría 
el 80% del almacenamiento total del país y esto equivale a 150 000 
hm? (Conagua, 2016). La Conagua, como autoridad responsable de la 
administración en materia de aguas nacionales y sus bienes públicos 
inherentes, ve la necesidad de conservar y, en su caso, mejorar la 
seguridad de las presas por medio de instrumentos regulatorios 
tendientes a mitigar los posibles efectos negativos o daños a terceros 
debidos a una descarga de sus obras de desfogue, o a su falla parcial 
o total, que pudiera provocar pérdida de vidas humanas o daño a 
zonas urbanas, infraestructura, y medio ambiente. El artículo 100 de 
la Ley de Aguas Nacionales establece que la Comisión Nacional del 
Agua elaborará normas, o realizará las acciones necesarias para 
evitar que la construcción u operación de una obra hidráulica altere 
desfavorablemente las condiciones de una corriente, o ponga en 
peligro la vida de las personas y la seguridad de sus bienes, o la 
integridad de los ecosistemas (NMX-AA-175-SCFI-2015). 


De esta forma, se ha propuesto una serie de tres normas que permite 
evaluar la seguridad de presas en México; sin embargo, la parte de 
estudiar los fenómenos extremos que inciden en estas estructuras 
hidráulicas es poco estudiada. En su primera parte, la norma 
establece los requisitos que deben cumplir las presas en operación en 
el territorio nacional para determinar su grado de riesgo o potencial 
de daños en la zona de posible inundación aguas abajo debido a la 
operación, falta de mantenimiento o de cuidado de las presas por una 
descarga, o por falla parcial o total de la presa (NMX-AA-175-SCFI- 
2015). 


En la segunda parte se detallan las inspecciones de seguridad de 
presas que son indispensables, y se realizan para determinar las 
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condiciones estructurales u operacionales de una presa o de alguna 
de sus partes (Ribeiro-Corréa, Cavadias, Clément, 8 Rousselle, 
1995). El tipo de inspección realizado dependerá de su propósito. 
Esta norma mexicana establece las especificaciones técnicas que 
deben cumplir las presas en operación en el territorio nacional para 
determinar el plan de inspección de presas (NMX-AA-175/2-SCFI- 
2016). 


En la tercera parte se formula la necesidad de contar con un Plan de 
Acción ante Emergencias (PAE) para coordinar las acciones que se 
deben realizar, a fin de proporcionar una notificación o advertencia 
para una evacuación en caso de emergencia; evitar o reducir la 
pérdida de vidas humanas y daños materiales o ambientales en el 
entorno geográfico de la presa preliminarmente en riesgo alto, o cuya 
clasificación formal de consecuencias por su potencial de daños sea 
alta. Esta norma mexicana establece los requisitos que se deben 
cumplir en la elaboración e implantación de un PAE para una presa; 
para su aplicación en caso de una descarga extraordinaria; la 
inadecuada operación de las obras accesorias; la eventual falla de 
alguno de los componentes de la presa, o la invasión de los cauces, 
de conformidad con lo establecido en la Ley General de Protección 
Civil (NMX-AA-175/3-SCFI-2017). Estas normas mexicanas son 
aplicables a los concesionarios, asignatarios, administradores, y 
usuarios formales e informales de las presas en operación dentro del 
territorio nacional y corresponde a éstos observar su cumplimiento. 


Con esta normativa, es claro que México trabaja hacia la prevención 
de desastres y la seguridad de presas. Aldama y Ramírez (1999) 
definieron la solución probabilística a la estimación de hidrogramas de 
avenidas de diseño, emplearon un análisis de frecuencias conjunto de 
gastos y volúmenes, ambos máximos anuales. Desarrollando con ello 
una solución bivariada (Aldama, Ramírez, Aparicio, Mejía, € Ortega, 
2006), para diseñar y revisar vertedores de 16 presas, cuyo periodo 
de retorno conjunto fuese de 10 000 años. En un estudio 
complementario, Campos-Aranda (2010) identificó los patrones de 
comportamiento estadístico en las soluciones bivariadas previamente 
referidas, desarrollando con ello un análisis crítico de las mismas. Lo 
anterior llevó a un procedimiento para diseñar o revisar vertedores de 
presas medianas o pequeñas ante la ausencia de datos hidrométricos 
del sitio, y contando con datos regionales de gastos máximos 
anuales. De aquí que en México, al ser tan vulnerable por fenómenos 
extremos, es necesario detallar estudios de seguridad de presas y la 
combinación de estudios multivariables. 
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En México, 224 ciclones tropicales impactaron en las costas entre 
1970 y 2015. En la temporada de huracanes 2015, el de mayor 
categoría que impactó fue Patricia (H4) (Conagua, 2016). 


Es importante mencionar que las condiciones y usos de suelo en las 
cuencas que se tenían en 2015, cuando ocurrió este huracán, no son 
las mismas que se tienen hoy día, por lo que alguna correlación 
obtenida con datos actuales y de 2015 podría no ser válida, y por lo 
tanto un pronóstico de volúmenes dentro de un embalse bajo estas 
condiciones sería poco confiable. 


Lo anterior lleva a plantear una hipótesis, donde la evolución de los 
volúmenes de almacenamiento en los mismos intervalos de tiempo en 
los que ocurre una lluvia huracanada podría describir el 
comportamiento hidrológico instantáneo de la cuenca de aportación. 
Se dice instantáneo, ya que las imágenes de los huracanes están 
disponibles en un intervalo de tiempo mínimo de 15 minutos. 
Asimismo, se acepta que el volumen de agua dentro del embalse de 
una presa tiene una correlación directa con la intensidad de 
precipitación y representa la respuesta física de la cuenca. Tal 
respuesta está condicionada por diversos factores, como cobertura 
vegetal, tipo y uso de suelo, y las características geomorfológicas del 
sitio, entre otros. Es decir, la precipitación es el origen del 
escurrimiento que finalmente ingresa a la presa; por esta razón, se 
acepta que el cambio de almacenamiento refleja las características 
fisiográficas, climatológicas y ambientales de la cuenca. 


Sin pretender proponer un modelo lluvia-escurrimiento, el objetivo 
del presente trabajo es correlacionar los valores de intensidad de 
precipitación durante un evento huracanado (calculados con un 
hidroestimador), con los volúmenes dentro de un embalse. Esta 
estimación del volumen de almacenamiento en presas durante la 
temporada de huracanes, empleando correlación canónica, puede ser 
una alternativa para revisar políticas de operación de presas en 
México. 


Metodología 


Análisis de correlación canónica 
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El análisis de correlación canónica fue desarrollado como una nueva 
técnica multivariada, en donde se utilizan los valores y vectores 
propios (eigenvectores) para analizar de modo directo gradientes 
multiescalas (Sun 8 Kim, 2016). Se trata de una generalización al 
coeficiente de correlación de momento de producto de Karl Pearson 
(Rice, 1972), que establece la proporción de variabilidad compartida 
midiendo el grado de covarianza entre distintas variables linealmente 
relacionadas (Pontasch, Smith, 8 Cairns, 1989; DelSole 8 Chang, 
2003). 


Mediante el análisis estadístico lineal de múltiples variables se 
establece la relación óptima entre dos conjuntos (Friederichs 8 
Hense, 2003): la correlación simultánea de p variables dependientes 
a razón de q variables independientes. Se lleva a cabo mediante la 
combinación lineal de ambos conjuntos, estableciendo vectores de 
peso que maximizan su coeficiente de correlación (Cavadias, Ouarda, 
Bobée, 8 Girard, 2001; Ntale, Gan, 8: Mwale, 2003). 


El primer conjunto agrupa las variables independientes X= 
(x1, X2, ..., Xp); mientras que el segundo conjunto agrupa las variables 
dependientes Y = (y,,y2,..., Ya). Ambos se correlacionan a través de 
las variables compuestas U y W, definidas como variables canónicas 
(Gabriel, 1972). Dichas variables canónicas son el resultado del 
producto de la transpuesta de los vectores de pesos a = (a,,a),..., Ap) Y 
b=(b,,b,,..,by) por el conjunto de variables dependientes e 
independientes, respectivamente, lográndose su combinación lineal, y 
que presentarán su máxima varianza al mostrar ortogonalidad entre 
ambas. 


O = (Gta mita) = 04%, + 09X7 += + ApXp = aTX (1) 


W = (b,, bz, ..., dp) | % | = by, + b2y2 +0" + bpyp = DYY (2) 


Conforme a lo anterior, se define la matriz de varianzas-covarianzas 
muestrales de tamaño p x q, que corresponde a ambos conjuntos de 
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variables (Ouarda, Girard, Cavadias, 8 Bobée, 2001), precisando que 
Cyx = Chy 


Y = das e) (3) 


Cvx Cry 


De tal forma que el coeficiente de correlación entre las variables 
canónicas se denomina coeficiente de correlación canónica y queda 
definido como: 


— LOU) GCANb 
corr(U,W) = Jvar(U)y/var(w)  JaTCxxa/bTCyyb 2) 


En la ecuación (4) se observan dos restricciones de normalización 
para determinar los vectores de peso que favorecen su maximización 
(Khalil, Ouarda, 8: St-Hilaire, 2011). 


var(U) = aTCyya = 1 (5) 


var(W) = b"Cyyb = 1 (6) 


La solución al problema de correlación canónica considerando las 
restricciones de normalización se reduce a maximizar: 


corr(U,W) = a" Cyyb (7) 


Dicha maximización se realiza empleando la función Lagrangiana de 
dos parámetros: 


L(a, b) = al Cyyb = ala Cia = 1) — Tb" Cryb — 1) (8) 


El proceso implica la obtención de los gradientes multivariados, 
partiendo de la ecuación (8) en términos de ambos vectores de pesos 
de solución: 
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OL(a,b) 
== = Cxyb = 2T,Cxxa =0 (9) 
222 = Cijya — 217Cyyb =0 (10) 


Ambos gradientes permiten establecer un sistema de ecuaciones que 
determina los vectores de peso a y b que maximizan la ecuación (7). 


Cxyb = 2T¡Cxxa (11) 


Copa = 27,Cyyb (12) 


Para la solución se consideran las ecuaciones (5) y (6) para 
multiplicar, respectivamente, las ecuaciones (11) y (12), utilizando 
los vectores transpuestos para su simplificación: 


(Cxyb = uba aa? = a Cb = 274 (13) 


(Ch,a = 2t,Cyyb)b" = pcia = 2ZT, (14) 


Desarrollando el segundo miembro de las ecuaciones (13) y (14), se 
obtiene el escalar que facilita la solución y garantiza las restricciones 
de normalización: 


(aTCyyb)y" = b"Chya = 1 (15) 


Ambas ecuaciones son equivalentes, definidas por el eigenvector 
A =21, =2r,; por lo tanto, se pueden reescribir las ecuaciones (11) y 
(12) de la forma: 


Cxyb = ACxxa (16) 


Capa = ACyyb (17) 
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Reduciéndose el sistema, al despejar el vector a de la ecuación (16) y 
sustituyéndolo en la ecuación (17), se tiene: 


CA Lose b| = ACyyb 
xv [7 CxxCxv YY 


Crv Cry CxxCxyb— Pb =0 (18) 


De manera análoga, al despejar el vector b de la ecuación (17) y 
sustituirlo en la ecuación (16) se obtiene: 


Cxx Cxy Cpy Cya— Pa =0 (19) 


La solución de las ecuaciones (18) y (19) permite obtener los 
vectores a y b, que representan los eigenvectores que maximizan la 
correlación de las variables canónicas para el conjunto de variables 
del problema. 


El análisis estadístico temporal y espacial de la información necesaria 
para llevar a cabo la correlación canónica se hizo conformando y 
ordenando las series temporales de variables, de tal forma que se 
generaran dos conjuntos, que correspondieron a las variables 
canónicas, y obteniéndose los eigenvalores y eigenvectores que 
permiten maximizar la correlación entre ambas, estableciéndose el 
grado de relación existente entre ambos grupos de variables. 


Huracán Carlotta 


El huracán Carlotta tuvo su inicio formativo la noche del 13 de junio 
de 2012, a partir de la depresión tropical núm. 3-E, a 560 km de 
Tapachula, Chiapas, avanzando hacia la costa de Oaxaca, alcanzando 
categoría II en las cercanías de Puerto Angel. Al desplazarse hacia el 
oeste-noroeste incursionó en la región montañosa del estado de 
Guerrero, degradándose en la mañana del 16 de junio a depresión 
tropical, con ubicación a 80 km al noreste de Acapulco, iniciando así 
su proceso de disipación (SMN, 2012). 
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Para establecer la influencia temporal y espacial del huracán Carlotta 
conforme a la trayectoria descrita en las cercanías de los estados de 
Guerrero y Michoacán (Figura 1) sobre las presas de almacenamiento 
en esta zona, se obtuvo la información de los boletines de 
seguimiento emitidos por la Administración Nacional Oceánica y 
Atmosférica (Pasch 8 Zelinsky, 2012). A partir de esta información, 
se construyeron series temporales de las variables de la posición 
geográfica que corresponden a la longitud (Long.) y latitud (Lat.); 
además, se utilizaron las características de magnitud del fenómeno, 
que corresponden a la presión barométrica (PB) y velocidad de sus 
vientos (VV) (Tabla 1). 


Figura 1. Trayectoria del huracán Carlotta del 13 a 16 de junio de 
2002 en los estados de Guerrero y Michoacán, México. 


Tabla 1. Registro de los boletines de seguimiento del huracán 
Carlotta por parte del Servicio Meteorológico Nacional, en 
coordinación con el Centro Nacional de Huracanes. 


E Lat. | Long. PB vv 
Boletín Fecha/hora 
si a mb km/h 

1 13-06-2012 / 22:00 | 9.4 | -92.4 |11005.33| 55 
2 14-06-2012 / 04:00 110.3 | -93.3 |11003.67| 65 
3 14-06-2012 / 07:00 |10.6| -93.4 |1002.67| 70 
4 14-06-2012 / 10:00 |110.9| -93.6 |1001.67| 75 
5 14-06-2012 / 13:00 |11.3| -93.8 |1000.67| 80 
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6 14-06-2012 / 16:00 [11.8 | -94.2 | 999.67 85 
7 14-06-2012 / 16:00 12.1| -94.5 | 998.17 85 
8 14-06-2012 / 22:00 [12.5| -94.7 | 995.67 | 100 
9 15-06-2012 / 01:00 [12.8 | -95.0 | 993.83 | 110 
10 15-06-2012 / 04:00 [13.2 | -95.3 | 993.33 | 110 
11 15-06-2012 / 07:00 [13.5 | -95.6 | 992.17 | 110 
12 15-06-2012 / 10:00 | 14 | -96.0 | 989.67 | 130 
13 15-06-2012 / 13:00 [14.4 | -96.2 | 986.33 | 135 
14 15-06-2012 / 16:00 [14.8 | -96.3 | 981.33 | 165 
15 15-06-2012 / 19:00 [15.5| -96.6 | 976.33 | 165 
16 15-06-2012 / 22:00 [15.9 | -97.2 | 972.86 | 150 
17 16-06-2012 / 01:00 [16.4 | -97.9 | 977.00 | 130 
18 16-06-2012 / 04:00 [16.6 | -98.4 | 980.00 | 120 
19 16-06-2012 / 07:00 [17.2 | -99.3 | 985.33 75 
20 16-06-2012 / 10:00 [17.5 | -99.5 | 994.67 55 


La variación de los volúmenes en las presas de almacenamiento de la 
zona de impacto y disipación se hicieron evidentes por la influencia 
del huracán Carlotta (Figura 2, Tabla 2). La revisión del 
comportamiento histórico en el Sistema de Seguridad de Presas de la 
Comisión Nacional del Agua permitió la selección de las presas Ing. 
Carlos Ramírez Ulloa (El Caracol) y Valerio Trujano (Tepecoacuilco) 
en el estado de Guerrero, así como las presas Infiernillo y Tepuxtepec 
en Michoacán, para desarrollar el análisis de correlación canónica, 
utilizando además información proveniente de imágenes satelitales 
para el pronóstico de intensidad sobre su ubicación geográfica (Zhao, 
2015; Rana, Renwick, McGregor, 8 Singh, 2018). 
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Figura 2. Ubicación geográfica de las presas de estudio en los 
estados de Guerrero y Michoacán, México. 


Tabla 2. Registro de almacenamiento diario registrado durante el 
huracán Carlotta por el Sistema de Seguridad de Presas. 


Nombre oficial Ubicación Capacidad de almacenamiento diario (hm?) 
(nombre común) Long. Lat. 13 14 15 16 
Ing. Carlos Ramírez Ulloa | -99,995 | 17.948 1012.6 989.6 987.2 969.3 
Valerio Trujano -99.466 | 18.296 11.91 11.84 11.8 11.77 
Tepuxtepec -100.23 | 20.000 283.28 283.28 282.84 282.94 
Infiernillo -101.89 | 18.272 2 817.94 2 817.94 2 827.65 2 872.65 


Imágenes satelitales 


Se identificó el comportamiento de la evolución geométrica sobre 
ambos territorios estatales y de manera particular en la posición 
geográfica de las presas de estudio utilizando imágenes digitales 
GOES-13, las cuales fueron proporcionadas por el Departamento de 
Percepción Remota del Servicio Meteorológico Nacional (SMN) de la 
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Comisión Nacional del Agua (Conagua). Dichas imágenes son 
captadas en datum WGS84 sin tratamiento, con formato raw y 
comprimida en BZ2 por su tamaño; la nomenclatura asignada es 
acorde con la fecha y hora en que fue tomada, asignándole extensión 
PCX (Meza, Curiel, González, € Alvarez, 2014). 


Las imágenes presentan en la referencia geográfica, en su esquina 
superior izquierda, las coordenadas 36.4768* de latitud y -122.2590% 
de longitud. En su esquina inferior derecha son 14.1118* de latitud y 
-79.0817% de longitud, en una resolución de 8 bits, delimitada a 1 
280 pixeles para la longitud y 817 pixeles para la latitud. El tamaño 
de cobertura aproximada es de 2.36 km por lado. Respecto a las 
características del espectro infrarrojo (canal 4), presenta valores de 
brillo que van de O a 255, igual con su nivel digital (Meza et al., 
2014). La Figura 3 muestra la imagen infrarroja satelital de mayor 
intensidad del huracán Carlotta, que corresponde al boletín 9 emitido 
por la NOAA. 


do a 


Figura 3. Imagen infrarroja satelital de mayor intensidad del huracán 
Carlotta. 


Series temporales de información satelital 
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La lectura de los valores digitales de los pixeles de posición de las 
presas respecto de la imagen satelital y la obtención del pronóstico 
de intensidad de lluvia se llevó a cabo mediante el algoritmo de la 
técnica hidroestimador. Dicha técnica se desarrolló a partir de 
observaciones en la zona central de las grandes planicies de EUA de 
mayo a junio de 1995, a través de la estimación instantánea de lluvia 
en radares superficiales y las temperaturas del tope de la nube a 
través de imágenes satelitales, logrando establecerse una regresión 
estadística de dichas estimaciones (Scofield € Kuligowski, 2003). 


Ri 0410 ete (20) 


En la ecuación (20), la variable R representa la intensidad de lluvia en 
milímetros por hora; T es la temperatura del tope de las nubes en 
grados Kelvin 195%K < T <260%kK, que depende del brillo del pixel (B) 
en la imagen, considerando que si B > 176, entonces T = 418 — B; caso 
contrario, si B < 176, entonces T = 330 — (B/2). 


Se aplicó la ecuación (20) para los pixeles de las imágenes de satélite 
que corresponden a la hora y fecha de emisión de boletines, 
obteniéndose la serie temporal de la intensidad estimada en la 
posición geográfica de las presas consignadas (Tabla 2), así como 
identificar la cantidad de pixeles que conforman las siete zonas de 
evolución geométrica acorde con la intensidad de lluvia y que están 
asociadas con los rangos de la escala de niveles digitales (ND) de 
brillo mostrados en la Tabla 3, propuesto para identificar las zonas 
del huracán, denominándose (Z0) al centro del huracán y (Z6) la 
parte externa de la nubosidad (Molina-Aguilar € Gutiérrez-López, 
2017). La Figura 4 muestra la imagen satelital que corresponde al 
boletín de seguimiento número 9, interpretada en escala de 
intensidades, en la cual se observa la posición geográfica de las 
presas ubicadas en los estados de Guerrero y Michoacán. 


Tabla 3. Escala de nivel digital propuesta para la interpretación de 
las imágenes satelitales y la identificación de zonas de evolución 
geométrica durante el huracán Carlotta. 


Zona 
geométrica 


z6 Z5 z4 z3 z2 zi 


Brillo 
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(ND) 
. 0.00- 0.01- 0.08- 0.37- 8.66- 90.13- > 
Intensidad 0.01 0.07 | 0.36 8.65 90.12 107.81 | 107.82 
(mm/h) 
Y 4 
y) 


Figura 4. Imagen de satélite interpretada en escala de intensidades 
para las presas (O ) ubicadas en los estados de Guerrero y 
Michoacán, México. 


La identificación y lectura de las zonas de evolución del huracán con 
base en su extensión nubosa se hace en primera instancia definiendo 
un umbral para el nivel digital del brillo en 100. A partir de este valor 
se realiza el pronóstico de intensidad de lluvia medible por los 
sensores de las estaciones climatológicas convencionales O 
automáticas en superficie. Para relacionar dicha lectura con un sitio 
de interés, se efectúa la proyección de las coordenadas geográficas a 
la proyección Geos, en la cual se encuentra la imagen satelital, de tal 
manera que se establezca su posición terrestre con el pixel utilizando 
un script en lenguaje R conforme a las características del sensor del 
satélite ambiental geoestacionario operacional (GOES-13), que da 
cobertura a la república mexicana (Molina-Aguilar, Gutiérrez-López, € 
Trejo-Alonso, 2018a). 


Después, usando una rutina programada en el mismo lenguaje R, se 
transforma la imagen de su formato nativo con extensión .pcx a 
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imagen con formato .jpg, el cual se lee como matriz de pixeles, 
asignándole a cada uno de ellos un color conforme una escala de 
intensidades de precipitación, y generando un archivo de registro 
para aquellos pixeles que corresponden con la posición geográfica de 
estudio. La rutina, llamada get_info.R, realiza también el 
procesamiento estadístico (Molina-Aguilar, Gutiérrez-López, € Trejo- 
Alonso, 2018b), que define la cantidad de pixeles, así como los 
parámetros de tendencia central, dispersión y forma para cada rango 
propuesto (Tabla 3). Con esto se construyen las 13 series temporales 
de información asociadas con cada estructura hidráulica señaladas en 
la Tabla 2. Dos de estas series están asociadas con el volumen 
almacenado e intensidad de lluvia; cuatro caracterizan la trayectoria 
y magnitud del huracán; siete representan la evolución geométrica 
asociada con la intensidad de lluvia durante la duración del fenómeno 
medido en cantidad de pixeles. Se establecieron dos grupos de 
variables: el primero, considerando Vol y R; mientras que el segundo 
agrupó las 11 restantes (Long, Lat, PB, VV, Z6, Z5, Z4, Z3, Z2,Z1 y 
ZO). La Tabla 4 muestra la información recabada para la presa Valerio 
Trujano, utilizada en la aplicación del análisis de correlación canónica. 


Tabla 4. Series temporales empleadas para el análisis de correlación 
canónica en la presa Valerio Trujano. 


Boletín Vol. R Lat Long PB vv z6 Z5 z4 Z3 z2 zi ZO 
Hm? mm/h e a mb km/h Adimensional 

1 11.86 0.03 9.4 -92.4 1005.33 | 55 36372 | 50551 | 52942 | 61764 16075 418 1 

2 11.84 2.39 10.3 -93.3 1003.67 | 65 47826 | 53038 | 45315 | 45992 19149 1721 Z 

3 11.84 0.03 10.6 -93.4 1002.67 | 70 45185 | 54810 | 43988 | 45935 14860 201 0 

4 11.83 9.83 10.9 -93.6 1001.67 | 75 51452 | 49139 | 39521 | 41453 9425 566 7 

5 11.83 0.01 11.3 -93.8 1000.67 |80 45874 | 48817 | 36713 | 27730 9138 985 0 

6 11.82 0.00 11.8 -94.2 999.67 |85 39503 | 44514 | 31242 | 27372 7874 562 36 

7 11.82 0.00 12.1 -94.5 998.17 |85 36183 | 47761 | 37898 | 39057 22199 801 13 

8 11.81 0.00 12.5 -94.7 995.67 |100 39364 | 42100 | 38888 | 66248 50501 5258 937 

9 11.81 109.58 |12.8 -95 993.83 |110 45145 | 55709 | 55018 | 70499 50651 6504 1106 

10 11.80 0.14 13.2 -95.3 993.33 |110 59456 | 67322 | 58232 | 63325 34730 5635 552 

11 11.80 4.00 13.5 -95.6 992.17 |110 58740 | 70884 | 56574 | 58735 32042 2815 301 

12 11.80 1.42 14 -96 989.67 |130 51748 | 69517 | 59234 | 43809 34463 5638 262 

13 11.79 0.00 14.4 -96.2 986.33 |135 56299 | 63336 | 45354 | 37558 31407 6135 365 

14 11.79 0.00 14.8 -96.3 981.33 |165 42459 |51771 | 44943 | 45957 37088 3479 324 

15 11.78 0.00 15.5 -96.6 976.33 |165 41351 | 54034 | 44674 | 64601 44195 6136 1367 
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Para determinar los vectores de peso se utilizó la información de los 
primeros 15 boletines (Tabla 1), ubicando el valor de cada variable en 
la serie de tiempo con la fecha y hora correspondiente (Tabla 4). Para 
verificar la correlación establecida y definir la estimación se utilizaron 
los valores registrados en los últimos cinco boletines (Tabla 1), cuyas 
magnitudes se muestran en la Tabla 5. 


Tabla 5. Complemento de las series temporales de la presa Valerio 
Trujano para realizar la estimación del volumen de almacenamiento. 


Boletín R Lat Long PB vv Z6 Z5 z4 z3 z2 zl1 y 40) 
mm/h 2 mb km/h adimensional 
16 0.000 15.9 -97.2 972.86 150 41097 | 64428 | 67480 [83795 | 54153 | 11884 | 2850 
17 35.165 16.4 -97.9 977.00 130 60182 | 79405 | 75901 | 86634 |65371 | 13080 | 5253 
18 0.032 16.6 -98.4 980.00 120 62108 |93177 |69845 | 80492 | 44869 | 10929 | 8364 
19 0.020 17.2 -99.3 985.33 75 55208 | 91589 | 76647 | 78603 | 56924 19627 |2596 
20 0.000 17.5 -99.5 994.67 55 60035 | 79573 |67614 | 74739 [40626 |6054 |1504 
Resultados 


Una vez establecida la ubicación geográfica de las presas respecto de 
los pixeles en la proyección de la imagen satelital, se estimó, para 
cada una de ellas, la intensidad, empleando la técnica del 
hidroestimador en la fecha y hora de emisión de los boletines de 
seguimiento (Figura 5), poniendo atención en la influencia temporal y 


espacial para la precipitación en cada una de las zonas (Figura 5). 
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Figura 5. Series temporales de la intensidad de lluvia estimada en el 
pixel de posición para las presas de estudio al paso del Huracán 
Carlotta. 


Utilizando las series temporales de la Tabla 4 se planteó la 
conformación de las variables canónicas para la presa Valerio 
Trujano; U = a, Vol + a2R y V = b,; Lat + b2 Long + b3PB + ba VWV + 
bsZ6 + b6Z53+0b3Z4 + b3Z3 + bo Z2 + bio Z1 + bi: ZO, que 
permitieron definir la magnitud de los eigenvalores 21 = 0.999 y 22 = 
0.863, cuyos valores de correlación canónica son R; = 0.9991 y R2 = 
0.7455, respectivamente. 


Con el primer eigenvalor se resolvió la correlación, lo cual resultó en 
los vectores de peso solución a = (-1.000, -0.004) y b = (0.821, - 
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0.407029, 0.404, 0.109, 0.140, -0.087, 0.065, -0.159, 0.190, - 
0.230, 0.169), con lo que se pudo verificar el grado de 
correspondencia entre ambas variables canónicas, es decir, la 
correspondencia de las variables asociadas con el problema 
planteado. 


Una vez que se establece la correlación, se emplea la información de 
los boletines restantes (Tabla 5) para estimar el volumen de 
almacenamiento en función de las características de posición e 
intensidad del huracán, su cobertura nubosa en la región, así como la 
distancia implícita de la posición de la obra hidráulica respecto del 
centro del fenómeno meteorológico, tal como se muestra en la Figura 
6. 
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Figura 6. Correlación canónica de valores registrados (*”) para la presa Valerio 
Trujano, (a) y los valores de las variables canónicas para la estimación (O ) del 
volumen de almacenamiento (b). 


Para las presas Ing. Carlo Ramírez Ulloa, Infiernillo y Tepuxtepec se 
utilizaron los eigenvalores 41, = 0.996, A21 = 0.999 y %1 = 0.996, 
respectivamente, para la determinación de los vectores a y b de 
solución, que maximizan la correlación de las variables canónicas R?; 
esto, en cada uno de los casos congruente a las series temporales 
generadas para las variables descritas, lo cual permitió verificar su 
respectivo grado de correspondencia y llevar a cabo la estimación del 


Tecnología y ciencias del agua, ISSN 2007-2422, 10(6), 25-56. DOI: 10.24850/j-tyca-2019-06-02 
45 


2019, Instituto Mexicano de Tecnología del Agua 
Open Access bajo la licencia CC BY-NC-SA 4.0 


Tecnología y 


CienciaszAgua 


(https://creativecommons.org/licenses/by-nc-sa/4.0/) 


volumen de almacenamiento con la influencia del huracán Carlotta. 
La representación gráfica muestra el grado de correlación canónica 
existente en las presas de estudio restante (Figura 7), en función de 


los pesos que conforman cada una de 


respectivas (Tabla 6). 


las variables canónicas 
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Figura 7. Correlación canónica y valores de las variables canónicas 
para la estimación del volumen de almacenamiento: (a) presa Ing. 
Carlos Ramírez Ulloa, (b) presa Infiernillo y (c) Tepuxtepec. 


Tabla 6. Variables canónicas para las presas de estudio y el valor de 
la correlación existente entre las mismas. 
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Presa Correlación Variables canónicas 
CcanonicCa 
U = - 1.000 Vol - 0.004 R 
ai R?=0.9991  [|W = 0.821 Lat - 0.407 Long + 0.404 PB + 
PuJano 0.109 VV + 0.140 NZ6 - 0.087 NZ5 + 0.065 
NZ4 - 0.159 NZ3 + 0.190 NZ2 - 0.230 NZ1 + 
0.169 NZO 
U =- 1.051 Vol - 0.135 R 
E ie R? = 0.9917 W = - 2.157 Lat - 1.416 Long - 0.917 PB + 0. 
dad 621 VV + 0.096 NZ6 + 0.399 NZ5 - 0.516 NZ4 
- 0.182 NZ3 + 0.354 NZ2 - 0.333 NZ1 + 0.288 
NZO 
U = - 1.006 Vol - 0.023 R 
S z W = 1.157 Lat + 0.524 Long + 1.857 PB + 
Infiernillo R*= 0.9984 0.201 VV + 0.087 NZ6 - 0.071 NZ5 + 0.084 
NZ4 - 0.022 NZ3 + 0.090 NZ2 - 0.218 NZ1 + 
0.219 NZO 
U = - 1.038 Vol + 0.110R 
W = 4.973 Lat + 2.672 Long + 1.874 PB + 
A 0.315 VV - 0.057 NZ6 + 0.239 NZ5 - 0.024 
Tepuxtepes Dn NZ4 + 0.194 NZ3 + 0.235 NZ2 + 0.042 NZ1 - 
0.296 NZO 


Los valores de volumen estimado que será almacenado en el vaso de 
las presas, utilizando la correlación canónica, se calcularon a partir de 


las variables canónicas conformadas. 


Sustituyendo los valores 


restantes de las series temporales (Tabla 6) y que se concentran para 
su visualización y análisis (Tabla 7). 


Tabla 7. Comparación de los volúmenes registrado y estimado 
mediante la correlación canónica. 


ocn Valerio Trujano OS Infiernillo Tepuxtepec 
registrado estimado | registrado estimado registrado estimado registrado | estimado 
16 11.780 11.780 975.267 982.459 2857.650 | 2259.086 282.907 283.105 
17 11.776 11.740 973.029 986.017 2863.275 2766.194 282.919 283.165 
18 11.773 11.708 970.792 961.440 2868.900 | 2812.142 282.932 283.499 
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19 11.769 11.737 969.517 | 983.539 | 2870.977 | 2826.231 | 282.940 | 282.769 
20 11.767 11.726 970.168 | 992.655 | 2865.958 | 2810.828 | 282.940 | 282.393 


Las capacidades de almacenamiento al nivel de aguas máximo de 
operación son las siguientes: 39.386 hm? en Valerio Trujano; 1 
458.21 hm* en Ing. Carlos Ramírez Ulloa; 9 340 hm? en Infiernillo, y 
384.32 hm? en Tepuxtepec. Dichas capacidades serán el referente 
para verificar la diferencia del volumen estimado utilizando ACC 
respecto del volumen registrado (AVol), expresando en porcentaje 
(Y%AVol) la precisión en términos de dicha capacidad (Tabla 8). 


Tabla 8. Diferencia porcentual del volumen estimado respecto del 
volumen registrado acorde con la capacidad útil de cada presa. 


Valerio Trujano Ing: ir Infiernillo Tepuxtepec 
Boletín 

AVol % AVol AVol % AVol AVol % AVol AVol % AVol 
16 0.000 0.00% 7.192 0.49% |-598.564 | -6.41% 0.198 0.05% 
17 -0.036 -0.09% 12.988 0.89% -97.081 | -1.04% 0.246 0.06% 
18 -0.065 -0.17% -9.352 -0.64% | -56.758 | -0.61% 0.567 0.15% 
19 -0.032 -0.08% 14.022 0.96% -44.746 | -0.48% -0.171 -0.04% 
20 -0.041 -0.10% 22.487 1.54% -55.130 -0.59% -0.547 -0.14% 

Discusión 


El procesamiento de las imágenes de satélite en canal infrarrojo 
(Figura 3) —una vez georreferenciados los sitios de interés y aplicada 
la escala para la identificación de las zonas de evolución geométrica 
en el huracán (Tabla 3)— permite de manera automatizada la lectura 
simultánea de valores de nivel digital (Figura 4). El script en lenguaje 
R desarrollado ofrece la posibilidad de realizar las lecturas en la 
secuencia de imágenes que dan seguimiento al inicio, evolución y 
disipación de un fenómeno hidrometeorológico. Con la adecuación de 
los parámetros de los sensores en los satélites venideros la 
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metodología mantiene su vigencia, ofreciendo una posibilidad de 
análisis dinámico y ágil. 

Con dichas lecturas y auxiliados de la rutina complementaria, pueden 
construirse las series temporales de intensidad de lluvia empleando la 
técnica del hidroestimador. Dicha técnica requiere de la validación a 
las condiciones específicas de ciertas regiones de la república 
mexicana, situación que se está trabajando y se presentará en 
futuros trabajos, disminuyéndose con ello la incertidumbre en el 
pronóstico realizado. 


Como se aprecia en la Tabla 4, a pesar de la presencia del huracán en 
su etapa inicial desde el 13 de junio por la noche, la cobertura nubosa 
no ha alcanzado la posición geográfica de la presa Valerio Trujano. Es 
hasta 48 horas después que la aproximación del fenómeno a costas 
del estado de Guerrero y el posicionamiento de una región concreta 
de la cobertura nubosa, con las condiciones de temperatura y 
humedad, permite estimar precipitación en el sitio de 109.58 mm/h, 
lo cual muestra la severidad de la intensidad del huracán, 
disminuyendo después la precipitación hasta ser prácticamente nula. 
Ello, derivado de la dinámica atmosférica y del hecho de que la fuerza 
de los vientos desplaza la lluvia a otras regiones continentales de 
manera casi inmediata. En la Figura 8 se muestra la evolución en el 
número de pixeles (N) conforme se desarrolla el huracán, como se 
puede ver para las zonas Z3 y Z2, un gradiente importante de 
crecimiento hasta alcanzar un máximo relativo en la curva es 
precisamente el tiempo de máxima intensidad de lluvia estimada, 
como se consigna en el boletín 9. 
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Figura 8. Comportamiento del número de pixeles en las imágenes 
satelitales con fecha y hora de emisión de los boletines por parte de 
la NOAA. 


La zona izquierda (Figura 8) son los datos utilizados en la 
determinación de los vectores de pesos para establecer la correlación 
canónica (Tabla 6), mientras que la zona derecha con fondo gris son 
los valores empleados para desarrollar la estimación (Tabla 7). 


En la Figura 5 se observan las afectaciones producidas en distinto 
tiempo y con diferencia sustancial en la magnitud. No debe perderse 
de vista que la lectura sobre el pixel es representativa de una 
cobertura superficial de 2.36 km por lado, que son unos 5.6 km?, 


La información analizada hace evidente la influencia de la posición y 
evolución geométrica de los ciclones tropicales respecto del 
pronóstico de intensidad de precipitación (en función de la rotación de 
la nubosidad que acompaña el fenómeno hidrológico). Se observa 
que las presas ubicadas en latitudes mayores ven reflejado en 
primera instancia el efecto del huracán respecto de las de menor 
latitud. De igual manera, es evidente que espacialmente existe una 
alta variación de la precipitación en función de la latitud y elevación 
de cada presa. 


La consideración de mediciones en variables físicamente registradas 
como el volumen almacenado, en conjunto con las mediciones 
obtenidas por teledetección, permite una perspectiva más amplia de 
la interacción del fenómeno físico. Conforme se mejore la resolución 
de los equipos y sensores, la precisión en los satélites, las 
estimaciones, empleando imágenes de satélite, pueden mejorar de 
manera importante. No obstante, se puede apreciar que los valores 
producto de la medición satelital en la aplicación de la correlación 
canónica permite obtener coeficientes de correlación mayores a 0.99 
en todos los casos (Tabla 6). Lo anterior permite disminuir la 
incertidumbre en el error de valores estimados con esta técnica. 
Además, no se limita la posibilidad de incluir mediciones de variables 
meteorológicas que se encuentran ubicadas en las inmediaciones o al 
interior de las cuencas hidrológicas. 


En el caso de la presa Valerio Trujano, la máxima diferencia es de 
0.065 hm?*; para la presa Ing. Carlos Ramírez Ulloa, la máxima 
diferencia es de 22.487 hm*; en Infiernillo, 598.564 hm*; finalmente, 
en Tepuxtepec, 0.567 hm?, que respecto a sus respectivas 
capacidades de almacenamiento al nivel de aguas máximo de 
operación representan una diferencia de 0.16, 1.54, 6.41 y 0.15%, 
respectivamente, de dicha capacidad (Tabla 8). En términos 
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generales se puede apreciar que los resultados de los volúmenes 
estimados son menores que los volúmenes registrados y que sólo en 
el caso de la presa Ing. Carlos Ramírez Ulloa se da una tendencia a 
sobreestimar dichos valores. 


Se puede apreciar que en el caso de la presa Infiernillo se da la 
mayor diferencia en el pronóstico, siendo también que se trata de la 
mayor cuenca, con 109 443 km?, y también la de mayor elevación de 
cortina (151.50 m) de las cuatro, mientras que las presas con menor 
diferencia en el pronóstico corresponden también a las de menor 
capacidad y cortinas de mediana altura. 


Es importante mencionar que la correlación canónica establecida se 
deriva en exclusiva de los datos que coinciden con la fecha y hora de 
emisión de los boletines de seguimiento al huracán Carlotta, cuyo 
intervalo de tiempo son seis horas. No obstante, el SMN dispone de 
imágenes con resolución temporal de 15 minutos, lo cual ofrece la 
oportunidad de un análisis más detallado, siempre y cuando en el 
sitio de interés puedan obtenerse los volúmenes de almacenamiento 
con dicha resolución. 


La inclusión de valores paramétricos representativos de cada una de 
las cuencas en las presas de estudio y que son condicionantes de su 
respuesta hidrológica, como puede ser la geomorfología, tipo y uso 
de suelo, cobertura vegetal y la red de drenaje por citar algunos, 
incrementan la confiablidad de la correlación canónica obtenida con 
estas nuevas variables. Para mantener dicha confiabilidad es 
necesario validar los vectores de pesos periódicamente, producto de 
la variación progresiva de los valores paramétricos, en particular los 
asociados con los coeficientes de escurrimiento a consecuencia del 
cambio de uso de suelo, y la disminución o incremento de vegetación. 


Conclusiones 


La propuesta del análisis de correlación canónica para llevar a cabo la 
estimación de los volúmenes de almacenamiento en una presa sin 
utilizar características paramétricas de las cuencas, correlacionado 
con intensidades de precipitación, arroja buenos resultados. El 
pronóstico del volumen de almacenamiento afectado por la lluvia 
huracanada resulta con un bajo porcentaje de diferencia en términos 
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del almacenamiento al NAMO, lo cual permitiría tomar decisiones 
oportunas respecto al desfogue controlado para que el vaso tenga la 
capacidad de regulación. 


El análisis canónico como una técnica multivariada se puede 
implementar como una medida que incremente la seguridad de la 
infraestructura ante fenómenos extremos, y sobre todo los eventos 
anteriores o posteriores; se podría incluir la lluvia antecedente. 


Los resultados obtenidos establecen una oportunidad para explorar la 
inclusión de más variables en el estudio asociadas con la 
geomorfológica específica de las cuencas en las que se lleva a cabo la 
transformación  lluvia-escurrimiento. También podrían incluirse 
parámetros estadísticos del comportamiento geométrico del huracán 
como un fenómeno físico. Si se logra incrementar la precisión en la 
estimación, sería posible definir o revisar la política de operación de 
presas durante la temporada de huracanes en el país. 
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